Critical behavior of two-dimensional random hopping fermions with vr-flux 
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A two dimensional random hopping model with A''-species and vr-fiux is studied. The field theory 
at the band center is shown to be in the universality class of GL(4m,R)/0(4m) nonhnear sigma 
model. Vanishing beta function suggests delocalised states at the band center. Contrary to the 
similar universality class with broken time reversal symmetry, the present class is expected to have 
at least two fixed point. Large A'^-systems are shown to be in the weak-coupling fixed point, which 
is characterized by divergent density of state, while small A'^ systems may be in the strong-coupling 
fixed point. 



I. INTRODUCTION 



It is widely accepted that the metalHc states are un- 
stable in two dimensional (2D) disordered systems [Q. 
However, two typical exceptions have been known for a 
couple of decades. One is the system with strong spin- 
orbit coupling g] and the other is the integer quantum 
Hall (IQH) system [||-g|. Especially, IQH transition has 
attracted much interest, because, as shown by Pruisken et 
al. Q , what is responsible for the delocalization is a topo- 
logical term in the usual nonlinear sigma model (NLSM) 
description. Though his arguments clarified nonpertur- 
bative aspects of the delocalization phenomena in the 
IQH systems, it was not generally possible to calculate 
critical properties of the IQH transitions. In order to 
calculate them explicitly, a model of Dirac fermions with 
various types of disorder was studied extensively [p^pOJ . 
This model has also intimate relationship with a statisti- 
cal model [O-jlJ and with more physical systems such as 
t-J model [|14| or d-wave superconductor p5Ml7|. How- 



ever, many problems remain to be explored, since generic 
IQH transition is in the strong-coupling Hmit ^ , though 
some exact results can be obtained for models with rather 
limited realization of disorder, e.g. the model only with 
random vector potentials 0. 

On the other hand, Gade found a new universality class 
having a random critical point pq |. It is due to what is 
called particle-hole symmetry, which enlarges the symme- 
try of the universality class to the general linear group. 
The renormalization group equations were calculated for 
the NLSM specifying such universality class with p^ or 
without [ p^p9| time-reversal symmetry. Recently, it was 
recognized that the random flux model belongs to this 
universality class |20[| or to its supcrsymmetric version 

Recently, Hatsugai et al. (HWKM) |p2| proposed an 
interesting 2D model with the particle-hole symmetry. 
This model describes random hopping fermions on the 
square-lattice with 7r-flux. It should have close relation- 
ship with various kinds of systems mentioned above, since 
it is described by the Dirac fermions near the band cen- 
ter. Numerical calculations by HWKM strongly suggest 



a random critical point at the band center. This model 
is actually in sharp contrast to the similar model but 
with on-site disorder [g3[ which has only localized states. 
Especially, what is interesting is the power-law behav- 
ior of the density of state (DOS) such as p ^ E'^ with 
V dependent on the strength of the disorder. Their cal- 
culations show that the exponent v keeps positive even 
for rather large randomness. This behavior of the DOS 
seems to be in contrast to the pioneering work on the 
two-sublattice model by Oppermann and Wegner p4[ , 
where the two-particle Green function has singularity at 
the band center. Even in one dimension, this symmetry 
plays a role in the existence of the critical point [p5|-p9[ , 
which brings about a singularity of the DOS at the band 
center. Therefore, it is interesting to study the 2D model 
proposed by HWKM, especially paying attention to the 
behavior of DOS. 

In this paper, we study the 2D random hopping model 
with A^-species. The model with one species corresponds 
to the HWKM model. In order to take 7r-flux into ac- 
count effectively, we firstly derive the Dirac Hamiltonian 
and next apply the NLSM method developed by Gade 
|18[ . We show that the model can be described by the 
NLSM on the symmetric space GL(4to,R)/0(47tt,), which 
is just the class BDl by Zirnbauer PQl- We conclude that 
the band center is a random critical point. However, 
contrary to the class with broken time-reversal symme- 
try, the present model should have, at least, two fixed 
point, which are characterized by the scaling property of 
the DOS. It is shown that large A^ systems are in the 
weak-coupling fixed point, where the DOS diverges. A 
discontinuity is, therefore, expected at zero disorder. It 
is also conjectured that the A^ = 1 model proposed by 
HWKM, or more generally, small A^ systems are in the 
strong-coupling fixed point. 

In Section O, we introduce the model and derive its 



continuum limit. In Section HI, we calculate the gener- 



ating functional of the Green functions and discuss the 
symmetry of the model. In Section pA^ , we use the trick 
of auxiliary field and give the saddle point solution. Sec- 
tion M is devoted to the derivation of the nonlinear sigma 
model describing the Goldstone mode. Section n\ deals 
with the renormalization group. Summary is given in 



Section VII 



It should be noted that the following calculations are 
indebted to MM and especially to ^. 



II. MODEL 

We present a slightly generalized model of HWKM, in- 
troducing an additional internal degrees of freedom. The 
tight-binding Hamiltonian is defined by 



^-EE' 



\(y'^i'r,jrCjT + h.C. 



(2.1) 



(■'■j) 



where j = {jx,jy) and a = 1,2, •••,iV denote, re- 
spectively, a site on the square lattice and a species of 
fermions. The summation with respect to sites is over 
nearest-neighbor pairs. The hopping matrix elements are 
defined by 



t 



j+xa,JT 



i-YHS^r + 5t^,aT, 



^j-\-ya.JT ~~ ''^a 



dt 



y,(TT 1 



(2.2) 
(2.3) 



with X = (1,0) and y = (0,1). Here, Stx^ar and Sty^ar 
are assumed to be independent random variables. The 
iV = 1 model is just the one proposed by HWKM. As 
stressed in |22] , the point is that the model has random- 
ness in the hopping terms, which causes the particle-hole 
symmetry. 



H 



H 



for 



(-> 



jx+ja 



(2.4) 



Accordingly, if there is an eigenstate of energy e, there 
also exists a conjugate pair with energy — e. The zero 
energy states is, however, special, because there is no 
partner for it. 

Let us first consider the pure model without random- 
ness. In this case, particles with different species do not 
couple each other. In the momentum space, the Hamil- 
tonian for each species is 



H„^2t 



E/ t t \ ( cos{aky) cos{akx) \ ( Cak 
\^'yk^^ak+ko) \co<i[ak^) -cos{aky) 



— TT-^akj: <7r 
0<aky<7r 



Ccrfc+fco 



(2.5) 



where a is the lattice constant and ko = (0,7r/a). Eigenvalues are E = ±2tyJ coa^ akx + cos^ aky, and hence there 
exist two Fermi points given by kp = (±7r/2a, 7r/2a). In order to investigate the model near the band center, it is 
convenient to linearize the dispersion at the two Fermi points and take a continuum limit p3| , p2[ . This may be valid 
up to a certain energy scale specified by a cut-off A. The lattice operator takes the following form in the continuum 
limit, 



Cj^/a -- ipa{x) = V''^'^i)ai{x) + i^" ^"^'^2(2;) + i ^""*"-'''-0<T3(a;) + i ■'" ^^ipaiix), 



where x = aj. 



Now it is easy to derive the total continuum Hamilto- 
nian including effects of the random hopping. The result 
is 

H = (dPx-^'^'H-^, (2.7) 

,*jv) with *t ^ 



where *t ^ (*{,■• 
(V4i'V'l2''^L'V'l4)' and 






Ti = Tip + TLd , 



hp = ttf^id^, 



V — 2at, 



K'a 1 *"^ 



(2.8) 



(Wd),,, 

Here, Tip and Tid describe the pure and the disorder 
Hamiltonian, respectively, and /x = 1,2. V^^, and U^^^, 
{j ~ 1, ... ,4) are independent random variables associ- 
ated with the Fourier components of St^^aa' and 5ty_a-a' , 
assumed to obey the Gaussian distribution. 






T>[V\e 90 J ^ Tt,'' 

V[U\e 90 J ^ aa'l 



Matrices a's and 7's are defined by 

a\= ax® T3, a-i — (J-},® 1, 
a3 = l«)T3, a4 = cr3(g)ri, 
7i = 1 (g) T2 , 72 = ^2 (8) n , 

73 =0-1 (8) T2, 74 = cr2<X'l, 
7 = CTl (g) Ti . 



(2.6) 



(2.9) 



(2.10) 



The particle-hole transformation Eq. (2.4) is expressed 
by 7 as 7_ff7 = —H . In the case of the on-site disor- 
der p3| , disorder potentials in the continuum limit com- 
mute each other and hence they are trivially diagonalized, 
while in the random hopping case, four 7's are noncom- 
mutative. Moreover, there appear extra disorder terms 



described by matrices a's caused by multi-species. When 
N = I, a-terms disappear and this Hamihonian reduces 
to the one derived by HWKM. 

In order to apply the same techniques as Gade's to the 
present model, it may be convenient to switch into the 
chiral basis, where 7 is diagonal. To this end, make a 
transformation, 7 — > U^jU, where 



/ 1 



U = 



2V2 



\ 



i 1 + il — i 1 — i 

l+i -1-i 1-2 -1+i 

-I -i l + i 1-i -l + i 

\-l-i -1-i 1-i l-^y 



(2.11) 



The matrices defined in Eq. (2.10) are then converted to 



ai = CTS (g) T2, a2 = CTl (X) T2, 

a3 = l(g)T2, a4 = -cr2®ri, 
7i = -o'3®Ti, 72 = -o-i«)ri, 
73 = -l(g)ri, 74 = -0-2 <8)T2, 
7 = -1 (X1T3. 



In the new basis, the Hamiltonian is real and symmetric 
due to the time-reversal symmetry. Moreover, the Hamil- 
tonian has now only off-diagonal elements in r-space, as 
is expected for systems with the particle-hole symmetry. 
The T-space is th us as sociated with the ±-sublattices in 
the lattice model (b.ll). 



III. REPLICA METHODS 

In this section, we will firstly derive the connection 
of the Green functions between the lattice and the con- 
tinuum theory |23], and secondly develop the continuum 
theory using the replica method. 

A. Green function 

Single-particle (e.g. retarded) Green function of the 
lattice system is given by 



Gaf^{xa,x'<j';z)=-i / dte" {0\i}^a{x,t)i}l ^{x' ,0)\0) 



{xaa\{z-H)-^\x'a'/3). 



(3.3) 



In Eq. (3^), terms breaking translational invariance, 
i.e., off-diagonal Green functions Gap with a ^ fi are 
neglected. 

Since we are mainly interested in the DOS for the lat- 
tice system, 

p{E)^+- hm VlmGL(ja,j(T,-E;±zc^), (3.4) 



let us define the following combinations of the Green 
functions for the continuum theory. 



G {x) = N^ 2, Gaa {x(T, Xa, 



Zk) 



a CK— 1 



(2.12) K^^\x,x') = Y, E Gaa'{xa,x'a',zk) 



a.a' q.q' — 1 



y.Ga'a{x a ,xij,zki), (3.5) 



where fc specify the retarded and advanced Green func- 
tion. (See Eq. (3^) for more precise definition.) 



B. Generating functional of the Green functions 

Since the Hamiltonian is real and symmetric in the chi- 
ral basis, we can apply the real bosonic replica method. 
Generating functional of the Green functions is intro- 
duced as 



Z= fv^t 



-s-s^ 



where 



s = lJd'xJ2--^k'<^ai^k-n)<i>t 

a,k 



(3.6) 



(3.7) 



GL{j<J,fa';z) = -i dte'-(0|c,,(i)ct ,(0)|0), (3.1) with k^{p,q) 



where z = E + ius, Cja{t) — e^^'cjo-e"'^* and |0) denotes 
the fermion vacuum. This Green function is expressed in 
terms of the continuum fields as 

GL(j<T,jV;z)^e'^(---'+^-^') 
4 
X ^ e*^["°(=""^')+'"°(2^-^')lG'aa(a;cr, x'a'; z), (3.2) 

a=l 

where constant vectors in oscillating phases factors are 
n = (0, 1,0, 1) and m = (0,0, 1, 1), and where 



Zk = {~r+'E + i-r+'u; = {-y+^E ~ Skiu, 

Sk = -isgnlmzfc = (-)^i. (3.8) 

Indices a — 1, . . . , m of scaler field ^^ denotes the replica. 
In what follows, to simplify the notations, we sometimes 
use a = (a, k) and a' — (a', k'), and hence, e.g., $^ = $5. 
Fixing a, the field <I>a is multi-component with respect 
to the species '$^ ee (*$ia, . . . , *$^a, . . . , '$7Va) and 
to Dirac indices '<I>o-5 = {^aad) with a — 1, ... ,4. Ss 
is a source term which will be introduced momentarily. 
Green functions with quenched disorder is expressed by 



G''{x) = lim Sfe V V($<,„5(a;)$^„a(a;)), 

m—*{) ^ — ^ ^ — ^ 
(7 a— 1 

4 

i^^^'(x,x') = lim SkSk' 






m^O 



(T,(t' a,Q' — 1 



($„ 



'a(x')$^aa(a:;)*<Taa'(a;)*<T'Q'a'(x')), (3.9) 



where (• • •) is the expectation value with respect to 
5'. Since the two components of the r-space in Eq. 
( ^.12 ) play an important role in Gade's argument p8|, 
let us introduce two fields cf)^ and 0~ explicitly, each 
of which has Dirac indices a — 1,2 of the cr-space, 

i.e., ($^ia,$a2a,$a3a,$<T4a) = {Kl~a^ 't'tla^ <t''^lh^ <t''^2a^) ■ 

The two fields reflect the particle-hole symmetry in 
the original lattice Hamiltonian, and are referred to as 
r = ± fields. To be more concrete, let us write down the 
Lagrangian with respect to (j)^ fields. 



= £p + £b + ^d 


+ A, 




p + -Cb = 


2^^ 

a 


(7 




(V+a 


V;a) 


\ ho 


'ho 

Zk 



Ca 



(3.10) 



^^—'oj:-^^j:i: 



a,a' j = l 



w ,(V+ , V"~) 



7, 



% 

a = -E44e 



a a 
a' a 



'^ (7 a 
K'a 



(3.11) 






d 



/i+, 

aa' 

K, 



era' 



(3.12) 



where 

ho = TmI^M' 
71 = f3, 72 



0-1, 73 = 1, 



74 



-iO-2, 



(3.13) 



and £p, C\^, C^, and £s describe pure-, breaking-, 
disorder-, and source-term, respectively. We set w = 1 
without loss of generality. The parameters in the source 
terms are symmetric had,'{x) = ha'dix) and are, more 
exphcitly, given by h^'^'^,. 



Let us also divide the Green functions in Eq. (3^) into 
two kinds of contributions from cj) , and average them 
over disorder (|2.9|) 



g\x) = Y,gUx), 

r=± 

K {x,x') = ^ K„,{x,x'), 



(3.14) 



r.r' — lb 



where 



Gl{x)= lim5feEE(<^-a(^)</'ia(^)) 

(7 a— 1 

dZ 



lim 



"^0 5/iaa(^) 
K„,{x,x')= lim SfcS/c' y^ y^ 

r?i^O ^ — ^ ^ — ^ 



-fcfc' , 



h+=h-=0 
2 



cJiCT' a,Q' — 1 



('/'^'o'a(2^')'^la(2^)'/'^'aa'(2^)'^^'a'a'(^')> 



= lim - 



d^Z 



--0 4 dh^J^,{x)dhl,^{x') 



(3.15) 



/i+=/t-=0 



Overbars denote ensemble-average, defined by 



(O) = lim / V[V]V[U]V[<^]0 

xe-^-i?^'^S^vE,[(^i.^)^+(^L')^], (3.16) 



and hence, Z in Eq. ( J3.15 ) shou ld be considered as in- 
cluding the disorder term in Eq. (3.16). Integration over 
V converts C^ into the form of four-point interactions. 



^d = -:^E^fe^fe'E E 



i+ .(h- ,.(h+ _,0- ,.,. 



a, a' a.a' — l 



(3.17) 



As Gade discussed ||T8|, this Lagrangian has unique 
symmetry property, which is caused by the particle-hole 
symmetry of the original model. Namely, it is easily veri- 
fied from Eqs. ( |3.10| ), ( |3.12| ) and ( ^.17D that when Zk = 0, 
the Lagrangian is invariant under the transformation 



^a-a7 



'S^9 ^s(i>„^, 



(3.18) 



where g = g^'^, G GL(4m,R), and s = s^'^, is a matrix 
defined by s = Saa'SkSkk'- Since the symmetry-breaking 
terms are written as 






aar^aaj i 



iE 



J2 E ( V+.A0+, + V;aA0-J , (3.19) 



(T a — 1 



with A = Saa'i—y^'^^kk', GL(4to,R) is broken by £„ 
up to 0(4to). For non-zero energy, Ce implies that 
the symmetry group is 0{2m,2m), which is broken to 
0(2?ti)xO(2?7i) by the Cui- This suggests that all states 
with nonzero energies are localized. 



IV. AUXILIARY MATRIX FIELD THEORY 

A. Auxiliary fields 

In order to apply the well-known trick of auxiliary 
fields, let us first introduce two kinds of fields, 



where h ^ h^ + h~ , h ~ h'^ ~ h~ , and tr means the 
trace in terms of a and k indices, i.e., trO = J^a k ^aa- 
Let us further define 



X 



Pfa' == 2 ^ ^ {'tia-aKca' + ' 

(T Q— 1 

2 



50 

N 
50 

N 



s-^ ps^ + h + ujs, 



^r .90 1 1 .7 

Y — — S20-S2 — ih. 



(4.3) 



G oe.aT (J (y_a ) "> 



aaar (jOLa 



(4.1) 



a a— 1 

By the use of these fields, we can express 



where X and Y are matrices with indices X^^, and Y^^j^, . 
We are then ready to introduce auxiliary fields. Noting 
the following identities. 



N 



C^ + C,+U + -— tr (x2 + r2 
250 





—UJ tl ss 

-tr(/is3 


^ps"!, 

1 - 1 1 


U = 


.90 ^ 


(s2ps2)^ + (s2crs2)^ 



N 



tr (h'^ -h^ -LU^ + 2ujhs) , 



) /'p[g]I?[P]e'^K-''^[(«-^)'+(^-^)1=l, (4.4) 



(4.2) 



where Q and P are real and symmetric matrices, we have the following expression for the generating functional 
depending on auxiliary fields Q and P as well as (p , 



Z = f V[(j,+]V[^-]V[Q]V[P]e- 1'^''=^^'^^'^-^\ 



N 



C^ -— [tr (Q2 + p2) - LutT s{Q + iP) - LutT s{Q - iP)] +Cp + CE-tT [Qs^ ps^ +Ps^as^ 



250 



N 
50 



tr h+ (Q - iP) + tr h- {Q + iP) - ujtT s{h+ + h-) -2tr/i+/i"l , 



(4.5) 



where we neglect w'^ term, which vanishes after the replica limit (m -^ 0). Integrating out the fields (j> , we end up 
with the following action. 



V[Q]V[P]e-^\ 

S^ = (fx [tr (g^ + P^) - Lutr s(Q + iP) - ujti siQ - iP)] + — Tr Ln C 

J ^9q -'2 

/N 
(fx— hi h+ (Q - iP) + tr h' (Q + iP) - Lutr s(h+ + h-) - 2tT h+h^] , 
50 

where Tr means the trace in x-space as well as tr and the trace in a-space, and C is defined by 









'EkSaa'Skk'-iQ'a^a'-^Paa:) 



(4.6) 



(4.7) 



with Ek — (— )''+^_E. The Green functions are given by 



rn^a go l\ ' 



KfAx,x')= lim f;^) {{Q-r'^P)ll{x){Q-r^P)Z{x'))-^{l-6rr')5{x-x'). (4.8) 

m^o \Zgo/ \ / Zgo 



B. Saddle point solution 

Most dominant contribution to the action, which is as- 
sumed to be independent of the coordinate, should satisfy 
the saddle point equations, 

— [{Q ± iP) - LOs] 

2 [Ek - (Q ± iP)] 



:ix,x), (4.9) 



[Ek - (Q + iP)] [Ek -{Q- iP)] + V- 

where relation h^ = V'^Sap is used, and the factor 2 in 
the numerator comes from the summation with respect 
to a = 1, 2. If cj = 0, it is expected to be of the form 



P = 0, Q= {EQk+ Skq)5aa'5kk'- 



Then, Eq. (4.9) reduces to one equation, 

1 ^ 2{Ek - Q) 



m 



{Ek - Q? + V- 



:{x,x). 



(4.10) 



(4.11) 



At zero energy this equation is easily solved as, setting 
Ek ~ Eok = 0, 



A(e" 



1) 2 ^ Ae "0 , 



(4.12) 



where A is a cut-off mentioned in Section O. Therefore, 
DOS is given by 



, 2iVA _^ 
p{E ^0) e 'o . 



(4.13) 



This equation shows that, for arbitrary coupling go, non- 
trivial solution exists, and the GL(47n,R) symmetry is 
spontaneously broken. 



V. DERIVATION OF NONLINEAR SIGMA 
MODEL 

In this section, we set i? = and concentrate on the 
band center. We will divide the fields into longitudinal 
and Goldstone modes and derive their transformation 



law under GL(4?ti,R) transformation (3. IS). This will 



turn out to be of importance when we derive a gauge- 
independent effective action of the Goldstone mode inte- 
grating out the longitudinal mode. 



A. Transformation properties 

So far we have shown that the Lagrangian is invariant 
under the global GL(4m,R) transformation (3.18) when 
uj = 0. Hence the action for the auxiliary fields (4.(;) 
should keep the same invariance. To study the trans- 
formation property of the auxiliary fields, let us con- 



sider the Lagrangian (4.5). Due to the coupling term 



between (jj^ and Q, P, the transformation ( [3.1^ ) for the 
quadratic terms, (/'+„0+„ -> g4>ta<i>ta^9 and (t>aa(t>aa ^ 
s ^ g~^ s4>~^4'~^sg~^ s induce the transformation of Q and 
P, 



Q + iP^ s-^ \f^ S-- {Q + iP)s^ g-^ s"^ , 
Q-iP~^ s^gs-^Q-iP)s-^*-gs^. (5.1) 

In the same way, the source terms should transform as 



-it -1 i, + i -1 -i 
s 2 '■g '^s^h^s^g s ^ , 

s^gs 2/i s ^ gs^ . 



(5.2) 



Note that the saddle point solution is invariant under the 
transformation o G 0(4?7i). Therefore, we can decompose 
the Q and P fields as follows 

Q + iP = s^^ ^T^^s^L+s^T^^s^^, 

Q-iP = siTs-iL-s-i^Tsi, (5.3) 

where the field T{x) G GL(4m, R)/0(4m) describe the 
Goldstone mode, and longitudinal fields are parameter- 
ized as 

L± (x) =L^ + sq, L± = Ql (x) ± iPl (x) , 
Ql^-{R-sRs), Pl^^{R + sRs), (5.4) 

with a real and symmetric matrix R. The constant imag- 
inary shift sq as well as the last two parameterizations 
ensure [|33|that the integration over </> fields converges in 



Eq. (4.6). Moreover, it is easily shown that the degrees 
of freedom are equivalent: The two real and symmetric 
matrices Q and P arc now converted to the two real and 
symmetric matrices T {^T — T is one of possible gauges) 
and R. The following identities are hold; 



sLnS 



-L^. 



(5.5) 



Next, examine the transformation of T and L fields. 
First, we have to consider the action of g eGL(4TO,R) on 
T. 



gT{x) ^T'{x)o{T{x),g), 



(5.6) 



where T'{x) G GL(4m,R)/0(4m) and o{T,g) G 0(4m). 
It should be noted that the field o{T,g) is a nonlinear 
function of T{x). Applying the transformation (pj) to 
( |3.4D and considering Eq. (5.6), we have the following 
transformation laws 



T^gTo-\T,g), 



L- 



s-^o{T,g)s--L 



+.^^-1 



{T,g)s- 



sio{T,g)s-iL-s-io-\T,g)si. 



(5.7) 



Under the change of variables ( p.3[ ), the action (4.5) yields, 
Z= I V[T]V[L]I[L]e^^% 



5a 



(fx 



N 
9a 



tr L+L^ -LutrCTTy'^sH+s^ + ujtr CTT)s-H- s-^ 



— TrLnC 
2 



tr h+s^Ts-^L-^s^^^Ts^ + tr h^ s^i ^T^^s^ L+s^T^^s^^ - ojti- s{h+ + h^) - 2tr h+h' 



where I[L\ is a measure and C is converted into, after the gauge-transformation 



C 



-L+ 



-L- 



with 



{T-^df,T± ^Tdf^^r 



(5.8) 



(5.9) 



(5.10) 



It is easily verified that the action is invariant under the local 0(4?ti) transformation as well as the global GL(4m,R) 
transformation [p4| , except for the symmetry-breaking terms. The transformation laws of v and a fields are 



ovuo ^ + odaO \ 



(5.11) 



which means that w is a gauge field associated with the hidden local 0(4™) symmetry Q]. Actually, the hidden 
local 0(4m) symmetry reflects the fact that the parameterization of the symmetric space GL(4to,R)/0(4to) is not 
unique. Another comment is that the breaking terms are not g-invariant, but are gauge-invariant. Therefore, the 
total Lagrangian including the breaking terms is gauge- invariant, as it should be. 

B. Derivative expansion 

Integration over the massive mode yields an effective action for the Goldstone mode. To this end, let us divide the 
action into two parts. 



Z= j V[T\D[L]e-^-^'^^-''^^^'^\ 



Sl{L) 



(fx tri+i" 

2.90 



N 



TrLnA-Ln/[L], 



6S{T,L) ^ — TrLn(l + A-iS)+ [(fx— \~ujti- CTTyhH+s^ + ujtr CTT) 3-^-3- ^ 
2 J 2go I 



d^x^ 
90 



tr ft.+s5Ts"^L"s"^*Ts5 + tr h^ s^^ ^T-^3^ L+ s^T-^s^^ - ujtr s{h+ + h^) - 2tr h+h" 



where 



A = 






B = 



-s 2 7^(wp - a^)s2 



S2 7^(u^ -|-a^)s 2 



(5.12) 



(5.13) 



(5.14) 



It is reasonable to define effective Lagrangian for T field as St = {SS{T, L)), where (• • •) means the expectation value 
with respect to Sl- Since we are now interested in long distance behavior of the Goldstone mode, we will expand 
these nonlocal Lagrangian up to quadratic order in fields, and will derive the effective NLSM for the Goldstone mode 
up to 0(iV°) for large N system. 

1. Longitudinal mode 



In this subsection, we derive the leading Lagrangian of 
the longitudinal mode. In order to get an 0(A^") NLSM 



of the Goldstone mode, it is enough to derive the La- 
grangian of the longitudinal mode up to 0{N). Hence, 
the measure term in Eq. (5.12) can be neglected. The 
first term in Eq. (5.12) is expanded as 

tr L+L" = tvL^L^ + qtr s{L^ + L^) -4mq'^, (5.15) 

where the last term can be neglected, since it vanishes 
after the replica limit. The second term, TrLn A, can be 
expanded as 



order with respect to the derivatives. In this leading or- 
der, V does not appear because it is not gauge covariant: 
It enters the fourth order in the form of covariant deriva- 
tive dfj_ — v^. Thus, in the second order, we can simply 
neglect v. 




Tr Ln A = Tr Ln Aq + Tr Ln (1 

1. 



An'Lo) 



(a) 



TrLnAo - Tr A^^Lq ~ -Tr A^^LqA^^Lo, 



(5.16) 



where 



A = Ao- Lo, 



An- 






Ln = 



L+ 

in 



(5.17) 



The linear term with respect to Lq in Eq. ( p. 16 ) can- 
cel out that in Eq. ( 5.15 ), as exp ected . Collecting the 
quadratic terms in Eqs. (5.15) and (5.16), we finally have 




Sl{L) 
where 



1 
Ad 



(fx{trd^,L^dij^LQ 



hLtrL^L, 



-^0 , 



N 
Girq'^ 



Uq^ 



(5.18) 



(5.19) 



(c) 

FIG. 1. Diagrams contributing (a) to order A'^ and (b), 
(c) to order A*'". Wavy line denotes a^, while thin and thick 
lines denote, respectively, the propagator of Aq and of L field, 
summarized in Appendix [A^. 



Using 



The derivation of this action is outlined in Appendix H. 



TrLn(l + yl"^S) 



Ti A^^B 



1 



TtA-^ba-^b. 



2. Goldstone mode 



as well as 
A-' ^ A-' + A-'LoA-' + A^^LoA^'LoA^\ 
it turns out that diagrams summarized in Fig. yj contribute to 0{N'^). Accordingly, we reach 



In order to derive the NLSM for t he Go ldstone mode, 
let us expand the action SS in Eq. ( 5.13| ) up to second 



(5.20) 



(5.21) 



5t= Jd^xl^tral - hr'a^ + ^c^gtr [CTT)-' + (*rr)]| 



90 



qivh+s^{T^T)s^ - qtr h- s^^T^Ty^s^^ ~ ojtr s{h+ + h-) - 2tr h+h' 



(5.22) 



where coupling constants are given by 



1 _ N db-4(l + 4m) 
1) 47r 47r 

- = ^ + 0{N-y 
c 47r 

db - 0.4593, dc ~ 0.1453. 



+ 0(iV-i), 



(5.23) 



The contribution from the diagrams (b) and (c) is denoted by d\y and d^ , respectively. An outline of the calculation 
is summarized in Appendix H. Some comments may be in order. Firstly, what is characteristic in this Lagrangian 
is that there appears tr ^-term as well as the principal term. As we shall see later, this term plays an important role 
in the scaling property of DOS. Secondly, at the leading order of the large N expansion, the coupling constant 1/c 



is 0, and in the next leading order, it is a finite positive values. This sign is quite relevant to the renormalization 
group flow, as we shall see momentarily. Lastly, this Lagrangian is manifestly gauge- invariant. If we fix the gauge as 
*T = r, the action becomes 



d^x 



Ab Ac 2go 



d^x- 



N r 
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gtr/i+s5T2s2 -gtrft,"s"5T-2s-2 - utr sQi^ + h^) -2tr h+h' 



(5.24) 



This is an action of GL(4m,R)/0(4m) NLSM, similar to the one for the two-sublattice model with broken time-reversal 
symmetry derived by Gade. The Green function is 



■^k 



N 



g1(x)= lim-gsfc[((T±^)^S>"^] 



k'' 



:{xy)=\va^^(^^^' (^[s'^T^s'^y'^^^^^^^^ 



(5.25) 



VI. RENORMALIZATION GROUP 



By the use of the relation |33] between the noncompact symmetric space GL(4r7T,,R)/0(4r?T,) and the compact 
symmetric space U(4m)/SO(4m), it is easy to write down the beta functions of the renormalization group |19|j36|] 



1 /I 



4 V 2 



Pb{b, c) = -nh^ - 7 o" + " P + o o" + 7" + '^ P' 



19 

7 

48 



/3c(6,c) 



1 
16 

i Pb{h,c) 
nb"^ 



1 /3 



43 



^-I^C(3)).^+(7 + .C(3) 



9 



—n 
2 



6^ - 0(5^), 



(6.1) 



where n = Am and C is the Riemann zeta function. Here 
we changed the sign of b of beta functions in Ref. ^91, 
because the present symmetric space is noncompact. We 
also took the normalization of b and c into account by 
calculating the renormalization constants explicitly up 
to two-loop order. The scaling equations are, therefore, 
given by, in the replica limit (m — > 0) , 



db 

dl " 


= 0, 


dc 
dZ" 


-r' 



0(6*) 



, (6.2) 



where I = \nL with L being the length of the system and 
coupling constants are rescaled as 6, c ^ ^, ^, whose 
initial values are given by 



bo 



Co 



N-{dh-dc) iV- 0.314' 
2 

dh 



(6.3) 



Let us now discuss the flow of the coupling constants. 
First of all, vanishing beta function for b implies that 



E — state is delocalized |18 . For E j^ 0, the universal- 
ity class 0(2?Ti, 2m)/0(2m)xO(2m) suggests that such 
states are localized, as discussed in Section [II. There- 



fore, only the band center is just on the critical point for 
the present model. Actually, this is in good agreement 
with the numerical calculation of HWKM for the N = 1 
model. 

To clarify the flow of c, let us first consider the case 
where N is large enough to ignore the 0(A^^"'^) in Eq. 
(^.23| ). In this case, b is quite small while c is positive. 
The scaling equation ( |6.2D then tells us that c scales to 
0. Therefore, large N systems behave similar to the class 
with broken time reversal symmetry ||l8| . Namely, since 
the C function is given by |19] 



C=^^b- 



_^ _ 1 (4 - 3n2 - n') b' 
c — no 64 



b b^ 3 ,, , 
2+7-16^ (™-0), 



(6.4) 



the DOS at the zero energy diverges under the change of 
the length scale; 



p (X q exp 



(dl 



(6.5) 



On the one hand, this may be likely, since Opper- 
mann and Wegner already derived a singular behavior 
of the Green function of the two-sublattice model such 
as l/N{d—2), using the large N expansion. On the other 
hand, the present divergence of the DOS suggests a dis- 
continuity with respect to 170 at go — 0, since DOS is 
exactly zero at the band center when qq = 0. 

However, the numerical calculation for the A^ = 1 sys- 
tem by HWKM suggests the convergence of the DOS 
even for rather large go. Therefore, small N systems, at 
least N = 1 system, should belong to a different fixed 
point. Actually in Eq. ( |6.2| ), there is a nontrivial zero 
for /3c at 6 = 6c ~ 3.087. Accordingly, if bo < be, c flows 
to 0, while if 60 > be, c flows to infinity, provided that 



Co > 0. The scaling equation (3.2) up to four loop order 
shows, therefore, the existence of a strong-coupling fixed 
point in addition to the weak-coupling fixed point c = 
mentioned above. Since c diverges at a certain length 
scale I = Ic, the DOS in the strong-coupling limit may 
scale as 



p a q exp 




ex q. 



(6.6) 



Namely, DOS is expected to be convergent in the strong- 
coupling phase. The initial value of the present calcula- 
tion is 60 = 2.92 for A^ = 1, which is still in the weak- 
coupling regime, to be sure, but lies quite near to the 
zero of the beta function. 

Of course, precise studies of this strong-coupling 
regime are beyond our scope, since we are base d on the 
perturbations: The initial values bo in Eq. (6.3) is valid 
only for large N cases. The zero of the beta function 
depends, moreover, on the order of the loop expansion. 
Namely, 6c is computed as 6c = 2 and 3.087 for two- and 
four-loop order, respectively, but no 6c for one- and three- 
loop order. However, numerical calculations by HWKM 
for the A^ = 1 system suggests the convergence of the 
DOS, which tells that such model cannot be in the weak- 
coupling fixed point. 

Therefore, we conjecture the existence of a critical N 
which separate the strong- and the weak-coupling phases. 
It is quite interesting to observe the weak-coupling one 
which is suggested in this paper for the first time, and 



to determine the phase diagram more precisely in the 
go-N plane by nonperturbative methods or by numerical 
calculations. 



VII. SUMMARY 

We have investigated the 2D random hopping fermion 
model proposed by HWKM, using field theoretical treat- 
ments developed by Gade. Starting from the tight- 
binding Hamiltonian, we have firstly derived the con- 
tinuum Dirac Hamiltonian which effectively describes 
the band center of the lattice model, and next con- 
structed the generating functional of the Green functions 
averaged-over the disorder. It has been shown that the 
symmetry group is GL(47tt.,R), which is spontaneously 
broken to 0(4m). Integrating out the massive modes, we 
have derived a nonlinear sigma model on the symmetric 
space GL(4TO,R)/0(4m.), describing the Goldstonc mode. 
The beta functions of the renormalization group show 
that the band center is a random critical point, where the 
density of state diverges for large N systems. This corre- 
sponds to a weak-coupling fixed point, which may share 
basic properties with the two sublattice model studied by 
Oppermann and Wegner. 

However, due to a nontrivial zero of the beta func- 
tion, it is likely that the small A^ system is in the strong- 
coupling limit. This fixed point is still critical but pre- 
sumably with convergent density of state. This fixed 
point may occur at a large value of the coupling con- 
stant, so that it is beyond our perturbative theory. It 
is, therefore, quite interesting to study the phase dia- 
gram more precisely by nonperturbative methods or by 
numerical calculations. 
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APPENDIX A: CALCULATIONS OF THE LEADING ORDER ACTIONS FOR THE LONGITUDINAL 

AND THE GOLDSTONE MODES 



The basic propagator Aq ^ is 



^oL'aa'(2;,a;') = ^aa'4fc' 



where, 






(Al) 
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d^P 



Jp-(x-x') 



(27r)2 p2 
tpp 



Px Py 

(27r)2 p"^ + q^ \ Py -px 



Jp-{x-x') 



(A2) 
(A3) 



Let us first consider the longitudinal mode in Section V B 1 , Main calculation lies in the quadratic term with respect 
to Lq in Eq. (5.16|), calculated as follows: 



fd^x fd'x' [2A2\x,x')A^\x',x)tysL+{x')sL+{x) - Aj^^,ix,x')A^^^,Jx',x)tyL,{x')L+{x)] +(+--) 
d^p f d'^p' 2{p-p' -q^)e-''^P-P'>'=' 



d^x' 



(2^)2 7(27r)2 (p2+g2)(p/2+^2) 



d^x [2tTL+{x)L^{x) ~ x'^tTdf,L+{x)d^,L^ix)] 



(A4) 



where we used Eq. ( p.5[ ) as well as a local field approximation Lq {x') ^ Lq {x) + {x' — x)^df_iLQ {x) + ^{x' — x)f_i{x' — 
x)iyd^d^LQ{x). The first tr term in the above equation as well as the first one in Eq. ( 5.15| ) give the mass term in 
Eq. (5.18), while the second tr term becomes the kinetic term 



In order to calculate the diagrams in Fig. |i| for the Goldstone mode in Section VB2 , we need the propagator of 
the longitudinal mode, given by 
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^ fl'ip)L'^wi-p)) = {s'l^Sl+C^i 



bj p2 



hL 



r± kk' 

-'O aa' 



'•ip)LVwi-p)) = -'^'^'' (-^aie'f' +'5„i'e'0 ^ 



(A5) 



The last equation follows from Eq. ( |5.5| ). By using these formulas, we have the following expressions contributing to 
the action (^.22) from each diagram in Fig. IT], 
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d a; (tr ( 
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(c) = 2iVd(l + 4m) 



d^p f d'^p 



g2(2p-p'+p2) 



(27r)2 J (2^)2 (p2 + g2)3(p/2 + ^2)[(p _ p/)2 + /,^] 



d xtr a,, 



(A6) 



Coefficients of these terms (integrals over p) are in principle dependent on go- To illustrate this, let us consider the 
first one, for example: Since it is convergent, we can integrate it over the whole 2D momentum space, giving the 
value N/AiT. However, for the present model, the cut-off is a physical parameter, denoting the scale under which the 
linearization procedure of the dispersion relation is valid. Accordingly, if we introduce the cut-off to the (convergent) 

integral, it turns out to be dependent on go through the saddle point solution q such as ^ ( 1 — ^^2, 2 ) = ^(1"^ ^)- 
However, for quite small go, the exponential in this expression can be neglected again. 
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